Exchange and spin-fluctuation superconducting pairing in cuprates 
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We propose a microscopical theory of superconductivity in Cu02 layer within the effective two-band 
Hubbard model in the strong correlation limit. By applying a projection technique for the matrix 
Green function in terms of the Hubbard operators, the Dyson equation is derived. It is proved 
that in the mean-field approximation d-wave superconducting pairing mediated by the conventional 
exchange interaction occurs. Allowing for the self-energy corrections due to kinematic interaction, 
a spin-fluctuation d-wave pairing is also obtained. T c dependence on the hole concentration and k- 
dependence of the gap function are derived. The results show that the exchange interaction (which 
stems from the interband hopping) prevails over the kinematic interaction (which stems from the 
intraband hopping). 



I. INTRODUCTION 

Recent experiment on field-induced superconductivity 
in infinite layer cuprate compound CaCu02 has proved 
that superconducting pairing is confined to thft,copper- 
oxygen plane both for electron and hole doping. El There- 
fore to describe the mechanism of superconducting pair- 
ing in cuprates one should study interactions in one 
coppeiQOxygen plane. As was initially pointed out by An- 
dersono, strong electron correlations in copper 3d-states 
play the most important role in explaining antiferromag- 
netic (AFM) and insulating properties of the undoped 
plane. These correlations in the electronic liquid of doped 
electrons or holes could be also responsible for supercon- 
ducting pairing in the plane.u 

Electron correlations are usually considered within the 
framework of the Hubbard model or the so-called t-J 
model which can be derived from the Hubbard model in 
the limit of strong correlations!] Superconducting pair- 
ing mediated by AFM spin-fluctuations in the weak cor- 
relation limit of the Hubbard model has been proved 
by using different approaches: by applying AFM spin- 
wave modelta, within thepjfluctuation exchange approxi- 
mation (FLEX) (see, e.gErH and the references therein), 
a renormalization group technique (seeLfu and the refer- 
ences therein). A number of phenomenological models 
like the nearly AFM Fermi liquida was also proposed. 
■ In the strontpgrrelation limit a lot of numerical work 
has been doneElo though the obtained results are still 
controversial. For instance, ar4£>bust <i-wave pairing 
was observed for the t-J modeljLd while absence of the 
long-range order was reported for the original Hubbard 
modeLEJ Concerning analytical calculations, usually a 
mean- field-approximation (MFA) both fop-th&.t- J model 
(see, e.g.jla) and for the Hubbard modeltJilZI were em- 
ployed. Self-energy corrections were considered only for 



the t-J model within an involved diagram techniqueMI 
or by using theJXlori-type projection technique for the 
Green functions^. A numerical solution of the Dyson 
equation irJ13 revealed a non Fermi-liquid behavior in the 
normal state, while the d-w&ve superconductivity medi- 
ated by the exchange and spin-fluctuation pairing was 
observed. To clear up the origin of the controversial nu- 
merical results reported for the t-J and Hubbard mod- 
els and to elucidate the pairing mechanism in cuprates 
a comprehensive study of the Hubbard, or more general 
p-d model beyond the MFA is very important. 

In the present paper, we consider a microscopical the- 
ory of superconductivity in Cu£)2 layer within the effec- 
tive two-band Hubbard modelE£l in the strong correlation 
limit. We have proved, by a direct calculation of anoma- 
lous correlation functions in MFA that the d-wave pairing 
in the model is mediated by the conventional AFM ex- 
change interaction as in the t-J model. The retardation 
effects in the exchange interaction, which originate from 
the interband hopping with a large excitation energy, are 
negligible that results in the pairing of all the electrons 
(holes) in a conduction subband and a high T c propor- 
tional to the Fermi energy. Allowing for the self-energy 
corrections beyond the MFA, we consider also the spin- 
fluctuation pairing induced by kinematic interaction in 
the second order. The latter, acting in a narrow energy 
shell of the order of a spin- fluctuation energy u g , results 
in a lower T c . 

The paper is organized as follows. In Sec. ||, the gen- 
eral formalism and the Dyson equation for the matrix GF 
in Nambu notati on a re derived for the two-band Hub- 
bard model. Sec. Ill is devoted to the discussion of the 



superconducting pairing within MFA. The self-energy is 
calculated in Sec. [iV] within the self-consistent Born ap- 
proximation. Numerical results and their discussion arc 
provided in Sec. Conclusions are presented in Sec. VI. 
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II. GENERAL FORMALISM 

A distinctive property of cuprates is the strong anti- 
ferromagnetic superexchange interaction for copper spins 
that reaches a record value of 1500 K for transition metal 
compounds. It is caused by the strong pda hybridization 
t p d — 1.5 eV for the 3<i-copper states and the 2p-oxygen 
states and the small splitting energy of their atomic levels 
A p d ~ 3 eV. At the same time strong Coulomb correla- 
tions for the 3d states of copper, Ud — 8 eV, considerably 
increase the energy of two-hole 3d states and the lowest 
energy-Level at hole doping appears to be a Zhang-Rice 
singletO. These features of the electronic spectrum in 
Cu02 plane can be described within the framework of 
the p-d model HamiltonianE3: 

H = J2i €d dfad™ + € P C t c i°} 
ia 

Operators df a and cf a describe the creation of one-hole 
d and p states at sites i of the square lattice in the 
Cu02 plane with the energies and e p = td + A p d, 
respectively. Because of the large Coulomb correlations, 
only singly occupied 3d states are taken into account: 
df = df a (l — nf _ a ). For the bonding oxygen orbitals 
the Wannier representation is used that results in the p-d 
hybridization parameters Vy — 2t p dVij with the coeffi- 
cientsEj v$ = Vjj ~ 0.96, v\ = Uj j± a , — —0.14, 



V 2 



■j J±Qx±a ~ —0.02 for the single-site hybridization Vq 
and the hybridization between the nearest v\ and the 
next nearest v 2 neighbor sites, respectively, where a x / y 
are the lattice constants. Since the single-site hybridiza- 
tion is quite large, Vo — A p d , and Vq » one has 
at first to diagonalize the single-site part of the Hamil- 
tonian ([!]) and then to apply a perturbation theory for 
intersite hopping—pAs— a result of this cell-perturbation 
theory (see, e.g.|!9£i±a), one gets an effective Hubbard 
model with the lower Hubbard subband occupied by one- 
hole Cu-cZ-like states and the upper Hubbard subband 
occupied by two-hole p-<i-like singlet states 
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H.c.)}, 



(2) 



where we introduced the Hubbard operators: A" m = 
\in)(im\ for the four states n,m = |0), \<r), |2) = | fj,), 
a = ±1/2, a = — a. The energy parameters are given 
by Ei = id — A 1 an d E 2 — 2E\ + A respectively, where 
Id is a reference (renormalized) energy of the ci-hole, /x is 
the chemical potential, and-A ~ A p d is the renormalized 
charge transfer energy (seec£l). Here and in what follows, 
the superscripts 1 and 2 refer to the one-hole and singlet 
subbands, respectively. The hopping parameters can be 
written as tff = K a p Vij . The coefficients K a p depend 



on the dimensionless parameter t p d/A p d- Assuming the 
realistic value A p d = 2t p d, we get for these parameterscj: 
Kn ~ -0.89, K 22 ~ -0.48, K X2 ~ 0.83. Therefore, the 
effective hopping parameters for the nearest-neighbors in 
the model (||) are given by t e ff ~ \K 22 \ 2v\ t p d — 0.14t p d 
that results in a narrow bandwidth W = 8i e // ~ t p d- 
Since in (||) the charge-transfer gap A plays the role of 
the Coulomb repulsion U in the conventional Hubbard 
model, the effective Hubbard model (||) corresponds to 
the strong correlation limit due to a large ratio of A to 
the band width: A/W ~ 2 . Therefore at half-filling the 
Hamiltonian (||) describes the Mott-Hubbard insulating 
state (see Fig. la in Ref.El). 

In the present paper we have assumed a simple version 
of the p-d modelE3 which has only two fitting param- 
eters, t p d and A p d, in Eq. (|l|). Starting from a more 
general model which takes into account finite values of 
the Coulomb repulsion on c£-sites Ud, on p-sites U p , and 
between the nearest neighbor p-d sites U p d, and p-p hy- 
bridization tpp one arrives to the same effective Hubbard 
model (||) but with renormalized hoppingnarameters t"^ 
and charge transfer energy A (see, e. g.,E3o). 

The Hubbard operators entering (||) obey the multipli- 
cation rules: A" m A 4 fe ' = 8 m ^X" 1 and the completeness 
relation 
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The latter rigorously preserves the constraint of no dou- 
ble occupation at each lattice site i by any quantum state 
\in). 

To discuss the quasi-particle (QP) spectrum and su- 
perconducting pairing within the model Hamiltonian (||) , 
we introduce the two-time anticommutator retarded 4x4 
matrix Green function (GF) in Zubarev notationEj 



G ija (t-t') = ((X ia (t)\xl (t')» 



2T 



s — icu (t, — t 



') 1^^) G a (q, cj), 



(4) 



for the four-component Nambu operators X ia 

X\ a = (Xf Xf Xf Xf), (5) 

and Xi a obtained from (^|) by Hermitian conjugation. 
The GF (j§ can thus be written as a 2 x 2 supermatrix of 
normal, Gijo-(uj), and anomalous, Fij a (uS), 2x2 matrix 
components: 



(G) 



where the superscript (T) denotes the operation of trans- 
position. 

To calculate the GF (Q) we use the equation of motion 
method. Differentiation with respect to time t of the 
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GF (Q) and use Qf-.the Fourier transform result in the 
following equationo 



ojG ija (tu) = 6 l3 x+((Z ia \X]J u 



(7) 



where Zi a = [Xi a , H] . Assuming that the system is in 
the paramagnetic state, {Xf a ) = (Xf a ) , we get for the 
matrix X = {{^i<T,X} a }): 



/ X2 X3 

Xi X3 

Xl X2 

\X3 o xi 



(8) 



For the diagonal terms of the matrix we get xi = 
(Xf 2 + X[-) = n/2, xi = (X™ + X~f) = 1- X 2 where 
n = 1 + 6 denotes the hole concentration. The anoma- 
lous correlation function X3 — (A 4 02 ) describes single-site 
pairing and for the d-wave symmetry it vanishes 



X3 = (Xf) = (c a c, T ) = 0. 
Here we used the identity 

y-02 _ y-(H y J.2 _ 



(9) 



(10) 



which results from the definitions of the Fermi annihi- 
lation operators: Ci a — Xf 17 + 2a X[ 2 and the multipli- 
cation rules for the Hubbard operators, X^X®' 7 = 0, 

xx 2 Xf 2 = o, xf l x} 2 = xf. 

The chemical potential /z is calculated from the equa- 
tion for the average number of holes, 



= {N i ) = j2(xn+2{xn 



(ii) 



Now by using the Mori-type projection technique we 
write the equation of motion for the operator Z ir7 in 
Eq. (Q) as a sum of a linear part and an irreducible part 

orthogonal to it, zj* , which originates from the inelastic 
QP scattering: 



— [Xi a ,H] — Eii a Xia- + Z, 



y(ir) 



(12) 



The orthogonality condition {{z\ l j\ Xt }) =0 provides 



the definition of the the frequency matrix: 

Mia = ({{X i<7 ,H],Xl}). 



(13) 
(14) 



The frequency matrix ( J13| ) defines the zero-order GF in 
the generalized MFA. In the (q, (^-representation, its ex- 
pression is given by 



where tq is the 4x4 unity matrix. 



X 



(15) 



Differentiation of the many-particle GF (fjj) with re- 
spect to the second time t' and use of the same projection 
procedure as in ( |l2| ) result in the Dyson equation for the 
GF (H|). In (q, ^-representation, the Dyson equation is 



(<7 CT (q,u;)) 1 = (GK^Lj)) l -E ff (q,w). 
The self-energy operator S CT (q, uj) is defined by 



(16) 



f a (q,w) = E„(q,w)+i: (T (q,w)G^(q,w)f (T (q,a;) ) (17) 
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denotes the 



where T a (q,u>) 

scattering matrix. From Eq. ( |17| ) it follows that the self- 
energy operator is given by the proper part of the scat- 
tering matrix that has no parts connected by the single- 
particle zero-order GF (| 



IP p i 

g: 



(ir) I y{ir)U\(prop) ~-l 



\z^))l—x 



(18) 



The equations (JTJ), (|TJ) 5 and ( plf ) provide an exact 
representation for the single-particle GF (f|). Its calcu- 
lation, however, requires the use of some approximations 
for the many-particle GF in the self-energy matrix Jl8| ) 
which describes the finite lifetime effects (inelastic scat- 
tering of electrons on spin and charge fluctuations). 



III. MEAN-FIELD APPROXIMATION 

A. Spectrum in the normal state 

Let us consider at first the QP spectrum in the general- 
ized MFA given by the zero-order GF (|l5|) . The latter is 
defined by the frequency matrix ( |l3| ) which can be readily 
calculated by writing down equations of motion for the 
Hubbard operators in Eq. (j^) and performing necessary 
commutations in the matrix Aij a , Eq. ( |I4| ) . As a result of 
these calculations, which are presented at greater length 
in the Appendix A, the matrix (^) can be written in the 
form 



(A*v) 



(T) 




(19) 



where Cjiju and Ay CT are 2x2 matrices. The diagonal 
parts of the matrix defines the normal state QP spec- 
trum: 



(E 1 + A) X2 + of af 

{af)* E lXl +af 

Kfl K? 1 



+ (i-Mv,i (A gj.^r], 



(20) 



while the off-diagonal matrices determine the supercon- 
ducting pairing: 
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bf bf 
-bl 1 bi 1 



+ (1 - Sij)Vij 



t22 t21 

ii(r ijcr 

_ J 21 T ll 

ij(7 ijcr 



(21) 



In Eq. the quantities a"^ determine energy shifts 

(renormalization of the chemical potential), while the co- 
efficients define renormalized hopping parameters. 
The site-independent anomalous correlation functions, 
b" 13 , and the site-dependent ones, Lf^ , in Eq. fl2l| ) de- 
termine superconducting pairing in MFA. Explicit values 
of the coefficients in E q. (20l) an d Eq. (|2l]) are given in 
the Appendix A, Eqs. (|A7H ]A1S| ), 

The QP spectrum in the normal state of the singlet- 
hole Hubbard model (§) within, the MFA has been stud- 
ied in detail by Plakida et al.B3. Therefore here we give 
only the results of the calculations which will be used 
further in discussion of superconducting pairing in the 
model (||) . The QP spectrum is described by the normal 
state components of the zero-order GF (jlF 







(22) 



-G°(-q,-u;)(T) 
where the 2x2 matrix GF reads0 

G°M,u;) = {cjfo-E rT (<l)r 1 
The q-representation of the energy matrix is given by 

EM - ( W W ) . (24) 



X2 

Xi 



(23) 



The energy spectra for the unhybridized singlet and one- 
hole excitations are defined by the equations 



W2(q) = E 1 + A + af/ X 2 + V 22 (q)/X2 
oj 1 (q)=E 1 +al 2 / X i +V^(q)/ X i, 

while the hybridization interaction is given by 



X2W' G 



The effective interaction in (E5|), (126]) has the form 



^(q) = |EW( k -l)' 



(25) 



(26) 



(27) 



where K£@(q) denotes the Fourier transform of K^f a in 



Eq. (|20|). Here we take into account only the nearest- 
and the next-nearest neighbor hopping by writing 



Kq) = 2y> ,e-^ 



3^0 



^i7(q)+8 i / 2 7 , (q), (28) 



where 



7(q) = 2( cos *r + cosq y ), V(q) = cos ^ cos ^ . (29) 



As was shown in RefR the off-diagonal matrix elements 
of the zero-order GF can be neglected since they give a 
small contribution of the order of t\ 2 /A to the density 
of states (see Fig. 2 in RefHil). Therefore, the zero-order 
GF (23) can be written in the diagonal form 



G°(q,w) 



X2 /[w-n 2 (q)] 

Xi/[w-Oi(q)] 



(30) 



where the hybridized spectra fi Q (q) for singlet (a = 2) 
and for one-hole (a = l) excitations are given by 



2 ,i(q) = i[w 2 (q) +w x (q)] ± ~{[w 2 (q) - ^i(q)] 2 



21tia121 1/2 



aw 21 w^ 2 } 



(31) 



In these equations and in what follows, the spin depen- 
dence of the QP energies is omitted. 

To obtain a closed solution for the zero-order GF, the 
correlation functions entering Eq. (pp[ ) are calculated self- 
consistently. The energy shifts a^^depending on single- 
particle correlation functions are readily calculated by us- 
ing the spectral representation of the GF (j30"l). However, 
the calculation of the functions K^ a depending on the 
dens ity-d ensity and spin-spin correlation function (see 
Eq. ( |Al2] )j-jjpquires the use of some approximations. 

In Refs.OliU, the Roth procedure was used which in- 
volved uncontrollable decoupling of the Hubbard oper- 
ators at the same lattice site, as e.g., (Xf 2 X 20 ) = 
(Xf a X[ 2 Xj°) . Here in the calculation of the effective 

hopping parameters K^ a we apply a simple approxima- 
tion by neglecting the charge-charge correlations NiNj 
at different lattice sites i ^ j (as in the Hubbard I ap- 
proximation) but fully take into account the spin-spin 
correlations: 



( X2 ) 2 + (SiSj) 



(32) 



For a spin-singlet state without long-range magnetic or- 
der, the GF ( p0|) and the one-hole spectrum ([n]) do not 
depend on spin. However, the short-range AFM spin- 
spin correlations are very important and, as calculations 
in Ref.Ej have shown, the spin-spin correlation functions 
(SjSj) in (§2|) bring significant contribution to the renor- 
malization of the dispersion relation (|3l|). For large 
spin-correlations at small doping values one finds a next- 
nearest neighbor dispersion. With doping, by decreasing 
the spin correlations, the dispersion changes to the near- 
est neighbor one which-js conventional to the overdoped 
case (see Fig. 1 in Ref.c3). The corresponding Fermi sur- 
face appears to be large even for a small doping as it is 
shown in Sec. Therefore, the present model repro- 
duces the most important features of the electronic spec- 
tra of cuprates: insulating state of the undoped Cu0 2 
plane and a strong reduction of the band width of single- 
particle excitations in the doped metallic state with the 
doping dependent dispersion. Concerning a non-Fermi 
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liquid behavior and a pseudogap formation in the un- 
derdoped region one has to study more accurately the 
self-energy corxeptions as has been done within the t-J 
model in Refs.EE] 



B. Superconducting pairing 



Let us consider the anomalous part Ay CT! Eq. (21), 
of the frequency matrix. The site in depen d ent a noma- 
lous correlation functions b^P , Eqs. ( A13 )-( [A15| ), van- 
ishes for the d x 2_ y 2 pairing. It can be easily proved 
by considering the q-dependent representation of the 
correlation functions, as e.g., J2 m ^i Vim { X i 2x m) = 
Y,qV((i){X° 2 X°l} . While the interaction V(q) is in- 
variant under the permutation of the components q x and 
q y (in the tetragonal phase) the anomalous correlation 
function having the d x 2_ y 2-wetve symmetry changes the 
sign that results in vanishing of the sum over q. 

To calculate the site dependent contributions , 
Eqs. (|A16|)-([A18|), we have to derive an estimate of the 
anomalous correlation function (X® 2 Nj). The approach 
followed in Refs.Ej'liH for its calculation employed the 
Roth procedure which decouples the Hubbard operators 
at the same lattice site by writing the time-dependent 
correlation function in the form: (cjj,(t)|cj| (t')Nj(t')) — 
(X° l (t)\X t l2 (t')N 3 (t')). In view of the identity satisfied 
by the Hubbard operators: X^ = X^X^ for any in- 
termediate state I7), the decoupling of the operators at 
the same lattice site is not unique and the refore unreli- 
able (a feature already noticed in Refs.HLH). 

Here we perform a direct calculation of the correlation 
function (Xf 2 Nj) without any decoupling by considering 
the equation for the corresponding commutator GF 



L ij (t-t') = ((X? 



'it) I Nj(t'))). 



After differentiation with respect to the time t and use 
of the Fourier transform, we get the equation: 

(lj - E 2 ) Lijiu) ~ 26ij {Xf) 
+ J2 2at\l{{{XfX^\N j )) u -{{XfX^\N j )) u }. (33) 

In the right hand side, only the interband hopping con- 
tribution (which mediates the exchange interaction pair- 
ing) has been retained since the contributions from the 
intraband hopping integrals give only a small renormal- 
ization of the large excitation energy \E 2 \ — A ^> |t^ Q | . 
The statistical average (Xf 2 Nj) at sites i ^ j can now 
be evaluated from the spectral representation theorem: 



f-ilm 



duj 



1 - exp(-w/T) 



£ 2 < 



2 
■tin 



1 



u> — E 2 + is 



({{XfX^Nj)) 

} 



uj-\-ie 



(34) 



We consider below a particular case of hole doping when 
the Fermi level stays in the singlet subband: fi ~ A, and 
the energy parameters, E% ~ E\ ~ —A . The contri- 
bution to the above integral coming from 5{u> — E 2 ) can 
be neglected since it is proportional to exp(— A/T) <C 1. 
The contribution from the one-hole subband Green func- 
tion in Eq. (|34|) can be estimated from its equation of 
motion as follows 



~lm{{XfX^\N 6 )), 



5 m] (XfX° j ")5(uj-2E 1 ), 



Therefore it gives also an exponentially small contribu- 
tion exp(— 2A/T) <C 1 ■ The only non vanishing contri- 
bution in Eq. (M) comes from the Green function of the 
singlet subband: 



(XfNj)^-- ]T 2atZ(XfX^Nj) 



(35) 



Here the retardation effect, the frequency dependence in 
the denominator l/(cu — E 2 ) , is ignored since the inter- 
band excitation energy is much larger than the QP exci- 
tation energy in the singlet subband: j^l — A 3> \t 22 \ . 
The exchange interaction is usually considered in the two- 
site approximation, which is obtained equating m = j in 
Eq. ([35]). In this approximation, we get after summation 
over the spin a: 

U 12 

(XfNj) = {cuaiNj) = -^-2a (X? 2 Xf 2 ). (36) 

where we have used the identity for the Hubbard opera- 
tors, XpNj = 2Xf 2 , and Eq. @. 

This finally allows us to write the expressions of the 
anomalous component A?? for the case of hole doping 
as follows 



(37) 



This result recovers the exchange interaction contribu- 
tion to the pairing, with an exchange energy parameter 
Jij = 4:(tj 2 ) 2 /A. The anomalous component A\^ a for 
the case of hole doping can be neglected since its con- 
tribution to the gap equation is extremely small, of the 
order of A^(q)/A. 

In the case of electron doping, on the contrary, we can 
neglect the anomalous correlation function for the singlet 
subband, A^ 2 CT , while an analogous calculation for the 
anomalous correlation function of the one-hole subband 
gives 



A 
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(38) 



We have therefore proved that in the MFA the inter- 
band anomalous correlation function is of the order of 
t]j/A and it is proportional to the statistical average 
of the conventional electron (hole) pairs at neighboring 
lattice sites in the Hubbard subband which intersects 
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the Fermi level. As a consequence, the anomalous con- 
tributions to the zero- order GF, Eq. @, originate in 
conventional anomalous pairs of quasi-particles and their 
pairing in MFA is mediated by the exchange inter actio n 
which has been studied in the t-J model (see, e.g.,t2rEj). 
In view of this conclusion, the MFA nonzero supercon- 
ducting pairing imparted in the frame of the conventional 
Hubbard modelEj Ej can be inferred to stem from the 
exchange interaction which equals to Jy = 4t 2 /U in this 
model. The exchange interaction vanishes in the limit 
U — > oo , a feature whieh-explains the disappearance of 
the pairing at large ^Jl3~tL3 

To summarize the study of the MFA, let us write down 
a closed system of equations for the particular case of 
hole doping with the chemical potential in the singlet 
subband. The 2x2 matrix GF for the singlet subband 
writes: 



X? 

xr° 



(39) 



By taking into account the normal part of the GF in the 
diagonal form, Eq. (|30|), and the anomalous part of the 
frequency matrix for the singlet subband, Eq. ([^), we 
can write the GF (k39j) in the form: 



(40) 



Gf^i^oj) = X 2 {cof - fi 2 (q)f 3 - 4>f (q)ri}" 



where tq, t\, T3 are the Pauli matrices. Here we intro- 
duced the gap function in MFA for the singlet subband 
</> 22 (q) = A 22 (q) 1x2 induced by the exchange interac- 
tion. From Eq. (f37|) and equation for the GF (^0|) we get 
a self-consistent BCS-type equation for the gap function 
in the MFA: 



^ 2 (q) = 



1 

N 



E^(k-q)f^tanh^M, (41) 
^ 1 ^2£ 2 (k) 2T ' V ' 

where J(q) = 4J7(q) is the Fourier component of 
the nearest neighbor exchange interaction and £ 2 (k) = 
[ft 2 (k) 2 + ^(k)! 2 ] 1 / 2 denotes the QP energy. Analo- 
gous equations can be obtained for the electron doped 
case, n < 1, with the chemical potential in the one- 
hole band by considering GF for the Hubbard operators 
X® a X*?® 

The obtained results for the gap equation in MFA for 
superconducting pairing mediated by AFM exchange in- 
teraction are identical to the MFA results for the t-J 
model (see, e.g.,E3£j). However, we have derived these 
results from the original two subband Hubbard model (g) 
which has allowed an understanding of the role of retar- 
dation effects in the exchange pairing. We have proved 
that they can be neglected due to the large interband ex- 
citation energy in comparison with the QP intraband ex- 
citation energy. We also obtained a more general, three- 
site representation for the anomalous correlation func- 
tion ( p5|) which can be used to study the many-particle 
exchange pairing. 



IV. SELF-ENERGY CORRECTIONS 



A. Self-consistent Born approximation 

To study the finite life-time effects and renormaliza- 
tion of the QP spectra caused by inelastic scattering one 
should calculate the self-energy (|l^) induced by many- 
particle excitations. A general representation for the self- 
energy in the (r, (^-representation reads 

y \ -*-■*■( Mij*b>) W 1 (A9\ 

S -H-X [ £t(x )(w) _Mg)(-,)J* - ^ 



where the 2x2 matrices M and $ denote the normal and 
anomalous contribution s to t he s elf-energy, respectively 
(see Appendix B, Eqs. @, ©). 

Here we consider the self-consistent Born approxima- 
tion (SCBA) (or the non-crossing approximation) which, 
is proved to be quite reliable in studies of the t-J model.E£l 
In SCBA, the propagation of the Fermi-like excitation, 
X\(t), and Bose-like exc itat ion, By{t), (spin or charge 
- see Appendix A, Eq. ( |A4j )) in the many-particle GF 
in ( [l2| ) are assumed to be uncorrelated. In the diagram 
technique SCBA is represented by a skeleton loop dia- 
gram without vertex corrections. Therefore the SCBA 
results from the decoupling of the corresponding oper- 
ators in the many-particle time-dependent correlation 
functions for lattice sites (1 ^ 1', 2 ^ 2') as follows: 



{B v {t)Xi{t)B v {if)X 2 {ff)) 
(X 1 (t)X 2 (t'))(B ll (t)B 2 ,(t')}. 



(43) 



Using the spectral theorem, the SCBA results in the fol- 
lowing decoupling relation for the many-particle GF: 



Hoc +00 



((B V X 1 \B 2 ,X 2 )) U ^^ J J 



dwidw2 

bJ — U)\ — UJ2 



x JV(wi,w 2 ) ImdX^X^ImdBv \B V % 



(44) 



where N(u) u u> 2 ) = (l/2)[tanh(wi/2T) + coth(w 2 /2T)]. 
Within this approximation we obtain a self-consistent 
system of equations for the self-energy (|42] ) and single- 
particle GF (g). However, to get a tractable problem, 
we simplify this system by using a diagonal approxima- 
tion for the GF as has been assumed for the zero-order 
GF (|l]). Consequently, keeping only the diagonal self- 
energy components, M" a and in Eq. (42), we can 
solve the 4x4 matrix Dyson equation ([16]) in the form 
of two independent 2x2 matrices for GF for the sin- 
glet and one-hole subbands. The diagonal approximation 
gives the lowest order contribution to the self-energy in 
the small hybridization parameter (^ 2 /A). 

By taking into account the zero-order 2x2 matrix GF 
for the singlet subband, Eq.(^), we derive the following 
Dyson equation for GF (|39|): 
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Gf (q,a;) = X2H - [fi 2 (q) + M 22 (q, w )]f 3 

-[^ 2 (q) + $f(q, W )]f 1 }- 1 . (45) 

Here the normal, M 22 (q, w) , and anomalous, $ 22 (q, w) , 
parts of the self-energy for the singlet subband in the 
SCBA approximation can be written in the form: 

+00 

M 22 (q, w) = 1 X] / dwxJCW Wl | k) q k ) 

K —00 

-ilm [if| 2 G 22 (k, Wl ) + ^(^(k, wi)] J , (46) 
$ 22 (q, W ) = -^]T /"dwiiC^^^l^q-k) 



— Im [^| 2 F 22 (k,o;i) - K^F^ik,^)] 



(47) 



The kernel of the integral equations for the self-energy is 
defined by the equation 



if( ± )( W , Wl |k,q-k) = t 2 \u(k)\ 2 



dtl>2 



N{ui,u 2 ) 

UJ — UJ\ — U>2 



-Imx^^q-k,^) 



, (48) 



where z^(k) is given by Eq. (28). The spectral density of 
spin-charge fluctuations is defined by the corresponding 
dynamical susceptibilities which are given by the com- 
mutator GF 

xit } q, <*>) = Xs(q, w) ± Xc(q, w) 

= ((S q |S_ q )) w ± i«fl\r q |<HV_ q » w . (49) 

Analogous expression for the GF in terms of the Hubbard 
operators A^ "", Af , and the corresponding self-energy 
components can be written for the one-hole subband in 
the case of electron doping (see Appendix B, Eqs. (B6)- 
(©)■ 



B. Weak coupling approximation 

Self-consistent numerical solution of the coupled sys- 
tem of equations for the GF ( f45| ) and the self- 
energy (|46|), (E3) is rather complicated (see, e.g. calcu- 
lations in Ref.tHf). To estimate the role of the kinematic 
interaction in superconducting pairing and its contribu- 
tion to superconducting T c beyond the MFA, we study 
below a simplified approach based on the weak coupling 
approximation (WCA). In WCA, it is assumed that the 
behavior of an electron liquid is dominated by the inter- 
actions around the Fermi level and therefore the interac- 
tion kernel ( f48| ) at frequencies (u>,u>i) close to the Fermi 
surface (FS) can be factorized in the form 



K (±) (w,wi|k,q- k) 
-itanh(^)AW(k,q-k) 



(50) 



for < lo s <C W where uj s is a characteristic pair- 

ing energy and W is the band width. In this approxi- 
mation, the effective interaction is defined by the static 
susceptibility 

A( ± )(k,q-k)=t 2 d |Kk)| 2 

+00 

x / ^[ilm X W(q-k,c 2) ] 

J W2 K 

-co 

/(k)| 2 Re X (±)(q-k,c 2 = 0)] . 



= t 



pd 



(51) 



The WCA is suitable for the band which crosses the FS. 
For the another band, which is far away from the FS at 
an energy of the order of the band gap, uj\ ~ A, an inte- 
gration over w 1 in Eqs. ([46]) and ( |47| ) is straightforward. 

For a hole doped system, n = 1 + S > 1, the chemical 
potential is in the singlet band, fi ~ A, and we can write 
the dispersion relations for the two bands in the normal 
state as follows 



fia(q) = n 2 (q) + x 2 " 1 M 22 (q,w = fia(q)) 
~ A-/i + e 2 (q) -e 2 (q) , 



fi!(q) = fii(q) + Xx 1 Ml\^u> = fi^q)) 
~ -M + ei(q) =i -A + ei(q), 



(52) 



(53) 



where the MFA energy £l a (q) is defined by Eq. (31). 
The zero of the energy in the singlet band is fixed at the 
Fermi wave- vector e 2 (q_F) — . 

Integration over u>\ in ( [47| ) gives the following equation 
for the anomalous component of the self-energy 



$ 22 ( q , w ~ o) = -k* 2 s 2 ,M) + Ki^M) 



(54) 



The sum <S 2 a (q) for the singlet band at the FS is given 

by 

5 2 ,.(q) = 1 V A<"> (k, q-k)|f| tanh ^1, (55) 



N 

k 



2£ 2 (k) 



2T 



where <I> 22 (k) = </> 22 (k) + $ 22 (q, 0) x 2 1 is an effective 
gap function. The integration over k is restricted here to 
an energy shell around the Fermi energy of the order of 
a characteristic energy uj s of the spin (charge) fluctua- 
tions defined by the susceptibility ([l9]) . The quasiparticle 
energy is given by 



&(q) = [0 2 (q)+ I # 22 (q) | 2 ] 



21 1/2 



(56) 



The one-hole band lies below the FS at an energy of 
the order A 3> W. Therefore, integration over W\ in 
Eq. @) for the anomalous GF F xl (k,wi) can be done 
by neglecting for the QP excitation energy in this band 
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the dispersion in Eq. (53) as well as the superconducting 
gap: Fffrwt) ~^(k)/K 2 -A 2 ). 
This results in the estimate 



Kk)| 2 



L^(k-q). (57) 



In this equation, the spin-charge static correlation func- 
tion iic } (q) = (S q S_ q ) - {(SNJiN^) resulted from 
the integration over u>2 in Eq. (Eq). A simple estimate 
shows that the sum <5>i i<T (q) gives a small contribution, 
of the order (t ef f / &f4>l l (k) ^ 1CT 2 <^(k) , and it can 
be neglected in Eq. (|54|). 

By taking into account the contribution due to the 
exchange interaction in MFA, Eq. (|d|), the equation for 
the effective gap in the singlet subband can be written 
as follows: 



(k) 



(k,q-k) 



k 

$22 



2£ 2 (k) 



tanh 



g 2 (k) 
2T 



(58) 



In this equation integration over k for the exchange in- 
teraction is performed without any restriction on the 
QP energy, while for the spin-fluctuation contribution, 
A<->(k,q-k), the energy of the pairing quasi-particles 
is confined to a narrow energy shell of the order of uj s 
close to the FS as pointed out in Eq. j50|). 

Similar considerations hold true for an electron doped 
system, n = 1 + 5 < 1 when fj, ~ in (^) . In that case, 
the significant WCA equation involves the gap $^(q) 
and critical T c defined by the corresponding equation 
analogous to Eq. fl5S|). 

Thus, we have proved that in the strong correlation 
limit of the effective Hubbard model (0), the equations of 
the gap functions of the two subbands: singlet and one- 
hole, can be solved independently since the self-energy 
corrections (both the normal and anomalous compo- 
nents) from the another subband give a vanishingly small 
contribution due to the large interband gap A 3> W. 



V. NUMERICAL RESULTS AND DISCUSSION 

We start with an analytical estimation of the super- 
conducting T c mediated by the exchange interaction and 
spin fluctuations in the gap equation (|5^) for a hole 
doped case. The charge-charge fluctuations contribution 
is small due to their large energy as compared to the 
spin-spin fluctuations and they are neglected. Therefore 
the effective interaction mediated by spin fluctuations in 
Eq. (pl|) can be written in the form 



A s (k,q-k)=i 2 Kk)| 2 Xs (q-k) 



(59) 



where Xs (q — k) is the static spin susceptibility. It can 
be evaluated from a dynamical spin susceptibility sug- 
gested in numerical studiesEa as follows: 



X.M « i(S,S_,) . TT ^ My (60) 

There are two fitting parameters in the model: a short- 
range AFM correlation length £ of the order of a few 
lattice spacing and a cut-off energy of spin-fluctuations 
lu s of the order of the exchange energy J. The spin 
susceptibility ( |60| ) has a maximum value Xo(£) a t the 
AFM wave vector Q = (tt, tt) when 1 + 7 (Q) = 0. Its 
value is fixed by the normalization condition for n — 
l + 5> 1: 

i^( Si s i )=^E( s < 1 s - q )=!( 1 -^' ( 6i ) 



which gives xo{0 =3(l-<J)/(4w a C(f)) where C(f) 
(V^O E q {l + C 2 [l + 7(q)]}- 1 - At large £ Xo (0 

e 2 /in£. 

For an analytical estimation of superconducting T c we 
consider the linearized Eq. (|58l) where for the gap we 
assume the <i x 2_ y 2-wave symmetry in the conventional 
form: 



$ 2i (q) = ^(cOBfc - cosq„) = ^7?(q) 



(62) 

Then integrating over q both sides of Eq. (|5§|) multiplied 
by ??(q) we get the following equation for T c : 



1 



A^ 



1 



20 2 (k) 



tanh 



n 2 (k) 

2T r 



[Jr,(k) 2 + A s (4 7 (k))Vk) 2 



(63) 



For the exchange interaction J(q — k) = 4J 7 (q — k) 
integration over q gives (4J/N) ^ q 7 (k — q) r?(q) = 
J7y(k) 2 . In q-integration of the static spin susceptibility 
in Eq. (55^) we have taken into account that it is posi- 
tive and shows a strong peak at the AFM wave-vector 
q k = Q = (tt, tt). Therefore the sum over q was eval- 
uated at this wave vector: (1/A) ^ q ?7(q)Xs(q — k) ~ 
-^(k)E q 'Xs(q') = -»?(k) 3(1 - 8)/4u>„. The effective 
coupling constant for spin-fluctuation pairing in Eq.(|6^) 
is given by A s = (K22 t pd 2v 1 ) 2 3(l - 5)/4w„ ~ t 2 eSf /uj s . 

Now we can perform integration over the QP energy 
0,2 (k) in Eq. (^) by introducing an effective density of 
electronic states (DOS) for two interactions: 



^( e ) = 77E^ k ) 2<5 ( £ -^(k)), (64) 

k 

^/( £ ) = Jj E ^ ( 4 ^( k )) 2 § ( £ ~ ^(k)) • (65) 



They are normalized to unity since (1/iV) ^ k ??(k) 2 = 1 
and (1/N) ^ k ^( k ) 2 (4 7 (k)) 2 = 1 . The resulting equa- 
tion for T r takes the form: 



W-fj, 



de e 

— tanh 

2e 2T C 



x [JN d {e)+d(uj s ~\e\)X s N sf (e)}, 



(66) 



s 



where 0(uo s — |e|) = 1 for QP energy |e| < ui s and equals 
zero otherwise. Here W is a renormalized band width 
for the singlet subband, Eq. (|52|). 

Let us compare superconducting pairing mediated 
by the exchange interaction and the spin- fluctuations. 
While for the exchange interaction - the first term in 
Eq. (|66|), the integration over energy is performed over 
the whole conduction band, in the spin-fluctuation con- 
tribution - the second term in Eq. (|6^), the integration 
over energy is restricted to a narrow energy shell close 
to the FS: uj s < E F . Moreover, the effective DOS @ 
for the exchange interaction includes large contribution 
from states near the van Hove singularity at (tt, 0) points, 
while for the spin-fluctuation interaction the correspond- 
ing DOS ( |65| ) is strongly suppressed since the kinematic 
interaction which is proportional to f(k) ~ 8^i7(k) van- 
ishes along the lines | k x | + 1 k y \ = tt . The effect of vanish- 
ing of spin-fluctuation pairing close to the AFM Brillouin. 
zone, I ky> | ~\~ | ky | = tt, was first mentioned by Schrieffero 
in his remarks concerning the model of nearly AFM Fermi 
liquid. □ 

These specific properties of interactions result in a 
lower T c mediated by spin-fluctuations in comparison 
with the exchange interaction though the coupling con- 
stant are comparable: J ~ 0.f5 eV, while A s ~ 
t 2 eff /uj s ^ 0.27 eV for t eff ~ 0.2 eV and uj s ~ 
0.15 eV. To demonstrate this we solve Eq. (66) by ap- 
plying the conventional logarithmic approximation since 
r c Cw s C(i. Then we get the following estimate for T c 
mediated by the exchange interaction 



rpe 



1.14 yj ^{W - ju)exp(-l/J N d (S)), 



(67) 



where N d (S) is an average DOS (Q) for a hole concen- 
tration 6. The superconducting T^ x peaks at an optimal 

doping for the chemical potential /i = W/2 and being 
proportional to a large electronic energy /i = Ep even 
for a weak coupling reaches a high value. For instance, by 
assuming the model parameters: Ep = W/2 ~ 0.35 eV, 
J ~ 0.15 eV, and N d {5) ~ 2 (eV • spin)" 1 , we get 
J N{§) ~ 0.3 and T c ~ 170 K. 

For the spin-fluctuation pairing we obtain the BCS-likc 
formula 



T*J ~ 1.14 Ws exp(-l/A s A s/ (0)), 



(68) 



where N sf (0) is DOS (65) at the FS. For a large FS, 



close to the lines | k x \ + \k y \ = tt , it is rather small as 
pointed out above that results in a weak effective cou- 
pling. Therefore Tj?f , being proportional to cu s <C /i 
and having weaker coupling should be lower than T^ x , 
Eq. (|67|). It should be stressed that in comparison with 
phenomenological models for spin-fluctuation pairing, as 
e.g. in Ref.u, the coupling constant X s N s f(0) in Eq.(|6"s|) 
is defined by the original microscopic parameters of the 
Hubbard model (Q). 

Similar estimations can be done for the electron doped 
case, n = 1 + S < 1, where /x ~ 0. In this case, the domi- 
nating contribution to T c comes from the integral S\ . For 



a conventional Hubbard model with only one hopping in- 
tegral tij we will have the electron-hole symmetry and a 
same T c . In the reduced p-d model (^), the band width 
for the singlet subband (proportional to tfj ) and the 
one-hole subband (proportional to tjj ) and the corre- 
sponding coupling constants are different that results in 
different T c (5) curves as observed in experiment .u 

A substantial proof of the analytical results given 
above is provided by a numerical study of the original 
gap equation (|5^) using a direct summation in k-space. 
The numerical solution of the gap equation ( |58| ) was per- 
formed under the condition of no double occupation at 
each lattice site i by a quantum state \in) in the upper, 
n = 2, or lower, n = a, Hubbard subbands: 



(at AT) = (AT xn = o, 



(69) 



which automatically follows from the Hubbard operator 
multiplication rules. The condition is satisfied by the or- 
der parameter with d x 2_ y 2 symmetry for a square lattice 
because the anomalous correlation functions (A£ 2 Af^) 
or (A^A *^) having the d x 2_ y 2-wave symmetry change 
the sign under the permutation of the components k x 
and k y that results in vanishing of the sum over k in the 
equation: (Af 2 Af 2 ) = (1/JV) £ k {X^X S _\) =0. 

The discretization of the Fredholm integral equa- 
tion ( |58| ) results in an eigenvalue problem for the gap 
parameter $ 22 (q) of the singlet subband. The temper- 
ature at which the largest eigenvalue equals one, while 
the corresponding eigenvector shows d x 2_ y 2 symmetry, 
will be taken for the critical temperature T c . We note 
that the eigenvalues of the integral equation are discrete 
(since the kernel is compact) and real (since the kernel 
can be made symmetric). 

The numerical results reported below have been de- 
rived under the following parameter values: A pc i = 
2t pd = 3 eV, t eff ~ K 22 2 Vl t pd ~ 0.14t pd ~ 0.2 eV. 
For the exchange interaction, we assumed the value 
J = QAt e ff, usually considered in the t-J model. In the 
model for the static spin susceptibility ( p0| ) the antifer- 
romagnetic correlation length was equated to the typical 
value £ = 3 which was kept independent of S and for the 
cut-off frequency we took uj s = 0.15 eV. 

Fig. |l| shows the dependence of the critical temper- 
ature T c (in t e ff units) on the doping parameter 6. 
The kinematic interaction alone results in a lower T c 
(the solid line in Fig. |l|) as compared to the exchange 
interaction (the dashed line). The highest T c (dotted 
line) is obtained when both interactions are included in 
Eq. (^8|). The maximum values of T c are quite high in 
all cases. They vary from T™ ax ~ 0.12i e// ~ 270 K 
at optimum doping S op t ~ 0.13 in the highest curve to 
T max _ 0Mt ef f ~ 90 K at S opt ^ 0.07 in the lowest 
curve. There is a reasonable agreement of these values 
with the crude analytical estimates (|67[), (|68[). It is im- 
portant to note that the WCA results overestimate T c 
because of the neglect of the inelastic scattering which, 
strongly suppresses T c as was proved for the t-J model.tS 
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FIG. 2. Wave-vector dependence of the gap function 
$ 22 (k) over the first quadrant of the BZ at optimum dop- 
ing (8 = 0.13) for temperatures: T — (a), 0.5T C (b), and 
0.9T C (c). The circles plot the Fermi surface. The +/— de- 
note gap signs inside octants. The numbers in graphs show 
the maximum isoline value. 



At smaller antiferromagnetic correlation lengths, the nu- 
merical solution suggested a strong decrease of the con- 
tribution of the kinematic interaction to T c . This feature 
stems from the decrease of the parameter xo(£) which 
controls the magnitude of the electron-electron coupling 
induced by the spin-fluctuations. The parabolic behav- 
ior of T c with 5, showing a maximum at optimum doping 
S op t — 0.13, is similar to the experimental data. This 
behavior originates from a strong dependence of the den- 
sity of states on the doping, with a marked peak at the 
optimum value 5ont , as shown in Fig. 2 in RefHd for the 
Hubbard model (H). At doping values 8 — ► 0, the super- 
conducting T c vanishes. However, at low doping a phase 
transition to AFM state should occur and therefore we 
do not show T c {8) for 5 — > . Moreover, at low doping a 
pseudogap formation should be taken into account that 
further suppresses T c . 

Further insight is obtained from the analysis of the 
k-dependence of the gap l> 22 (k) inside the BZ at differ- 
ent temperatures. In Fig. || a typical wave-vector de- 
pendence of the gap is shown in the first quadrant of 
the BZ, (0 < k Xl k y < 1, in n/a units), at optimum 
doping (S = 0.13) at three temperatures, T — (a), 
T = 0.5T C (b) and T = 0.9T C (c). The occurrence of a 
superconducting state with d x 2_ y 2 symmetry is observed, 
with a q-dependence pattern of the gap shown by iso- 
lines. This dependence is much more complicated than 
that depicted by the simple analytical model (|62|). 

Consideration of the kinematic interaction term alone 



in Eq.(58) and the subsequent numerical solution at the 



specific optimum doping value S = 0.07, resulted in gap 
maxima inside the BZ, while inclusion of the exchange 
interaction as well and numerical solution at the higher 
optimum doping value 5 — 0.13 resulted in the shift of 
the maxima towards the (ir, 0)-type points of the BZ, 
which are close to the FS, as shown in Fig. ||. The 
particular behavior of the spin-fluctuation pairing is ex- 
plained by vanishing of the kinematic interaction along 
the lines \k x \ + \k y \ = ir , close to the large FS, as dis- 
cussed above. These results show that the inclusion of the 
exchange term is essential for the achievement of exper- 
imentally observed in cuprates wave-vector dependence 
of the gap and an optimum doping values. An addi- 
tional experimental evidence for AFM exchange pairing 
mechanism in cuprates comes from the T c dependence on 
the lattice constants or an external pressure. While in 
most electron-phonon superconductors T c decreases with 
an increasing applied pressure, in cupeates T c (at opti- 
mum doning) increases with pressure.E!l As was shown 
recentlyo, the AFM exchange pairing model for T c ( |67| ) 
describes quite well the experimentally observed T c de- 
pendence on the lattice constants in the Cu02-plane in 
mercury cuprates. 

Therefore, we may conclude that the most important 
pairing interaction in the strong correlation limit of the 
Hubbard model (^) is the exchange interaction, while the 
spin-fluctuation coupling mediated by the kinematic in- 
teraction results only in a moderate enhancement of T c . 
The same result has been obtaiiiad in studies of the t-J 
model beyond the WCA in RefJlj. However, the super- 
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conducting exchange pairing in the t-J model is realized 
already in the simplest MFA for the nearest neighbor par- 
ticles, while in the Hubbard model a more complicated 
three-particle anomalous correlation functions, Eq. (35), 
are induced by a dynamical second order interband hop- 
ping responsible for the exchange interaction. This ob- 
servation pay explain why numerical simulations for the 
t-J modelHI usually predict a much stronger. pairing ten- 
dency than for the original Hubbard modelli 2 ] since in the 
latter more complicated dynamical processes should be 
accounted for. 



VI. CONCLUSIONS 

A microscopical theory of the superconducting pair- 
ing has been proposed within the reduced two-band p-d 
Hubbard model (|^) with two microscopical parameters 
only, the p-d hybridization parameter t p d and the charge- 
transfer gap A p d- By applying the projection technique 
for the two-time GF in terms of the Hubbard operators, 
we obtained the Dyson equation for the 4x4 matrix 



GF (16) with the corresponding matrix self-energy (18|). 

In the zero-order approximation for GF, Eq. (40[), 
we have proved that the <i-wave superconducting pair- 
ing of conventional hole (electron) pairs in one Hub- 
bard subband occurs that is mediated by the conven- 
tional exchange interaction Jy = 4(iy) 2 /A caused by 
the intcxhand hopping as in the t-J model. Contrary 
to Refs.E3 : tZl, the anomalous correlation function in MFA 
(X® 2 Nj) = (ciia^Nj) was c alcula ted without any decou- 
pling approximation (see Sec. HI B| ) . We have also proved 
that retardation effects in the exchange interaction are 
negligible that results in pairing of all electrons (holes) 
in the conduction band and a high-T c proportional to the 
Fermi energy, Eq. (|67|). 

Spin-fluctuation pairing was studied in the second or- 
der of the intraband hopping within the self-consistent 
Born approximation for the self-energy (f47|). The ob- 
served spin-fluctuation (Z-wave pairing in the Hubbard 
model is mediated by the kinematic interaction which 
vanishes inside the BZ along the lines \k x \ + \k y \ = it and 
produces pairing only in a narrow energy shell of the or- 
der of spin-fluctuation energy close to the FS: lo s <C Ep ■ 
These specific properties of the interaction results in a 
lower T c mediated by spin- fluctuations, Eq. (K38), in com- 
parison with the exchange interaction, Eq. ( |67| ) . Numer- 
ical solution of the gap equation (^8j) confirms the ana- 
lytical estimations for the superconducting T c and repro- 
duces experimentally observed T c {5) dependence with a 
maximum at the optimum doping 5 op t — 0.13, Fig. [I]. 
The cZ-wave symmetry of the gap function with a com- 
plicated k- dependence, Fig. ||, is also proved by the nu- 
merical solution. 

It is to be stressed that, within the Hubbard model 
in the limit of strong correlations, the electron-electron 
coupling induced by the exchange and kinematic interac- 



tions are caused by the non-fermionic commutation re- 
lations for the Hubbard operators, and therefore no ad- 
ditional fitting parameters are needed to accommodate 
that interactions. These mechanisms of superconduct- 
ing pairing are absent in_the fermionic models as was 
pointed out by Anderson.EJ However, they are different: 
while the AFM exchange pairing originates in the low- 
ering of kinetic energy of electron (hole) pairs caused by 
the interband hopping as suggested by Andersoro, the 
spin-fluctuation pairing is of the conventional nature: the 
pairing is due to attraction in the cZ-wave channel caused 
by the spin-fluctuation exchange. 

The weak point of the argument is the derivation of 
the reported results in the weak coupling approximation 
for the superconducting equations. To substantiate the 
present results, more rigorous self-consistent numerical 
solution of the strong coupling Dysonj-aquations in the 
(q, w)-space, as done for the t-J model,t2l should be elalt. 
orated. In this connection we can mention papersEj'Efl 
where superconducting pairing in the Hubbard model 
was investigated within the non-crossing approximation 
for the self-energy in the 4-cluster model in the dynamical 
mean field theory. 
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APPENDIX A: FREQUENCY MATRIX 

To calculate the frequency matrix (|l^) one needs to 
consider the equations of motion of the Hubbard opera- 
tors: 

Zf = [Xf, H] = (E\ + A)X° 2 

+ £ (^BZ l xf 2 -2atfB^,Xr') 

-^A° 2 (4 1 Ar° + 2a4 1 Ar), (Al) 
z oa = [xf,H]=E 1 X° i 3 

-Y^Xr{tfxr + 2at\fXf)^ (A2) 
Zf = -{Zf)\ Z? = -{Z«°)y. (A3) 
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Here, B^ a , are Bose-like operators describing the num- 
ber (charge) and spin fluctuations: 



B iaa> — { X i 2 + X i a )Sa'a + X[ a 5 a 'g 



B laa> - <W - B 



21 

zcxcx' ! 



(A4) 

(A5) 
(A6) 



-2cxK 21 ((2-A^)X° 2 ), 



(A16) 
(A17) 



I% a = -(K 22 (X™N 3 ) - Kn((2 - N^Xf 2 )). (A18) 



According to Eq. (|10[), the anomalous correlation func- 
tions describe the pairing at one lattice site but in differ- 
ent subbands: (X? 2 Nj) = (X° l X^Nj) = (a^Nj). 



where we have used the completeness relation m|), and 
the definition of the number operator Ni, Eq. ( Jll| ) and 
the spin operators Sf = Y,a <jX r^ S i = x f 9 ■ 
Performing necessary commutations of the operators 
Zf 2 and Z® s in the matrix A ij(7 , Eq. @, we get for 
the quantities a"' 3 in Eq. (EQ) the following equations 



<4 2 = £ 



V im {K 22 {X 2S X^) - K n (X^°X9-)) , (A7) 



2a £ V m K 12 ((Xf °X% ) - (X^fXf 2 )) . (A8) 

For the renormalized hopping parameters -K^f in the ma- 
trix ( p0| ) analogous calculations give the results 

Kfl^K^-K^Xrxf), (A9) 
ifu (X,7 + 1 - n) - K 22 (X? 2 Xf >, (A10) 



if 11 



(All) 



where stands for the static charge and spin correla- 
tion functions 



1 



X ^ = -(N l N J ) + (S l S J ). 



(A12) 



In Eq. (|2l|), the site independent anomalous correlation 
functions are given respectively by 



.22 



V m {K 22 {{XfX° m 2 )-{XfX%)) 



2aK 12 {{X? 2 X%) + {XfX%)) } , (A13) 
V im {K lx ({ X yX%)-{xrX%)) 



E 



b 2 } = 



2aK 12 {(X^X^) + {XfX%)) } , (A14) 

~v im {K 22 ((xrx^ + ixfxZ)) 



2oK 12 ({xrX%)-{xrX%))} 



(A15) 



The site-dependent anomalous correlation functions are 
convenient to write as follows 



APPENDIX B: SELF-ENERGY 

The starting point for the calculation of self-energy cor 
rections is the (r, ^-representation of the self-energy (18) 



(t f ) = r W } (*) 1 4 r)t (*')» brop) x- 1 - (bi) 



The operators Z ( ' r) are obtained from Eqs. (|Al|)-(|A^) 
according to definition, Eq. (12). Neglecting the terms 
which contain X® 2 operators, whose contribution to the 
low energy dynamics of the system is assumed to be small 
(since it involves large charge fluctuations), we get: 

uf^,x? 2 - 2oi$8B% a ,xy\ m 



y (ir) _ / + 22 XD 22 



7 (ir) 



I ^ uu il ulJ icrcr'^l 

WUnll vOa o,l 2xu l 2 v cr'2\ 



2atX6B% al X? 



(B3) 



Here, SB°f a , = B°f a , - (B°f al } are Bose-like operators 
describing charge and spin fluctuations as follows from 
Eqs. (0)-(|A6|) 



The many-particle GF in Eq. (Bl) associated to the 
irreducible operators, Eqs. (B2), (p33|), give the following 
2x2 matrices M and $ for the normal and anomalous 
contributions to the self-energy, respectively: 



Mij„{u) ■ 



M 22 » MfHu) 
Mf 2 a {u>) M#» 




(B4) 



(B5) 



In the diagonal approximations for the GF the self-energy 
matrix elements Mij a (u>) (B4) and 3>ij a (u)) (B5) within 
the SCBA (0) in the (q, ^-representation are given by 
the equations: 

+ OG 

/dwiJr (+) (w,wi|k,q-k) 

k 



12 



-Im 



. (B6) 



i-oo 

4 tT (q,w) = — ^ /dwiJT^^wilkjq-k) 



-ilm 

7T 



> 2 ( - ) F» 2 (k, Wl )-P 1 (_) i ; ; u (k,wi) h (B7) 



where 



p(±)_ 

— 



p(±)_ 
— 



A 22 



±2(7^21^22 

2aK 21 K 22 ±K^ 

K\ x ±2aK 21 K 11 
2aK 2 iK n ±Kf 1 



(B8) 
(B9) 



The diagonal matrix elements M| 2 (q, oj), $ 22 (q, Lj), for 
the singlet subband are given by Eqs. (Eg) , (|47|) , while 
for the one-hole subband, Mj: 1 (q, a;), $^(q, ap, can be 
easily obtained from the matrices (B6), (p7|). 
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